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We study spin effects of quantum wires formed in bilayer graphene by electrostatic confinement. 
With a proper choice of the confinement direction, we show that in the presence of magnetic field, 
spin-orbit interaction induced by curvature, and intervalley scattering, bound states emerge that 
are helical. The localization length of these helical states can be modulated by the gate voltage 
which enables the control of the tunnel coupling between two parallel wires. Allowing for proximity 
effect via an s-wave superconductor, we show that the helical modes give rise to Majorana fermions 
in bilayer graphene. 
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Introduction. Graphene and its derivatives [1-4], such 
as bilayer graphene (BLG) and carbon nanotubes (CNT), 
have attracted wide interest due to its peculiar band- 
structure with low energy excitations described by Dirac- 
like Hamiltonians. Moreover, these materials are usually 
placed on substrates, which allows high control of its ge- 
ometry, doping, and placement of metallic gates [5-9]. 
Topological insulators were predicted for graphene [10], 
but later it was found that the intrinsic spin-orbit inter- 
action (SOI) is too weak [11, 12]. For BLG, first-principle 
calculations also show weak SOI [13, 14]. In an other pro- 
posal, topologically confined bound states were predicted 
to occur in BLG where a gap and band inversion is en- 
forced by gates [15]. Quite remarkably, these states are 
localized in the region where the voltage changes sign, are 
independent of the edges of the sample, and propagate 
along the direction of the gates, thus forming effectively 
a quantum wire [15-17]. At any fixed energy, the spec- 
trum inside the gap is topologically equivalent to four 
Dirac cones, each cone consisting of a pair of states with 
opposite momenta. 

The spin degrees of freedom in such BLG wires, how- 
ever, have not been addressed yet. It is the goal of this 
work to include them and to show that they give rise 
to striking effects. In particular, we uncover a mecha- 
nism enabling helical modes propagating along the wires. 
In analogy to Rashba nanowires [18], topological insu- 
lators [19], and CNTs [20, 21], such modes provide the 
platform for a number of interesting effects such as spin- 
filtering and Majorana fermions [22]. Here, the SOI plays 
a critical role, and in order to substantially enhance it, 
we consider a BLG sheet with local curvature as shown 
in Fig. 1. Two pairs of top and bottom gates define the 
direction of the quantum wire which is chosen in such a 
way that it corresponds to a 'semi-CNT' of zigzag type. 
In this geometry, the energy levels of the mid-gap states 
cross in the center of the Brillouin zone. A magnetic 
field transverse to the wire in combination with interval- 
ley scattering leads to an opening of a gap, 2A^, between 
two Kramers partners at zero momentum, see Fig. 2. As 
a result, the number of Dirac cones changes from even 




FIG. 1. A bilayer graphene (BLG) sheet with a fold at x = 
along the z-axis is placed between two pairs of gates that are 
set to opposite polarities ±Vb/2, inducing the bulk gap. There 
are mid-gap bound states, localized in transverse x-direction 
around x = 0. At the same time, they freely propagate 
along the ^-direction, forming an effective quantum wire [15]. 
An externally applied magnetic field B = Be x / breaks time- 
reversal symmetry. The spin-orbit interaction /3 is induced by 
the curvature of the wire, which is characterized by the radius 
R. In the insets we show the BLG structure in momentum 
(left) and real (right) space for a chosen chirality = 0. The 
edges of the BLG sheet can be arbitrary. 



(four) to odd (three), and the wire becomes helical with 
opposite spins being transported into opposite directions. 
In the following we derive the spectrum and its character- 
istics analytically and confirm these results by indepen- 
dent numerics. We also address the physics of Majorana 
fermions which emerge when the wire is in proximity con- 
tact to an s-wave superconductor. 

Curved bilayer graphene with SOL We consider a 
gated curved bilayer graphene with a magnetic field B 
(along the x'-axis) applied perpendicular to the direction 
of the fold (along the z-axis), see Fig. 1. We begin with 
a description of the bilayer graphene in the framework 
of the tight-binding model [3, 4]. Each layer is a hon- 
eycomb lattice composed of two types of non-equivalent 
atoms A\ (A2) and B\ (B2) and defined by two lattice 
vectors ai and 3.2. We focus here on AB stacked bi- 
layer, in which two layers are coupled only via atoms A2 
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and B\ (see Fig. 1) with a hopping matrix element t± 
(t± « 0.34 eV). By analogy with CNTs [3], we introduce 
a chiral angle as the angle between ai and the x-axis. 

The low-enegy physics is determined by two valleys 
defined as K = — K' = (47r/3a)(cos 0, sin (9), where 
a = |ai|. The corresponding Hamiltonian in momentum 
space is written as 

H = hv F {k x a 1 +r 3 k z a 2 ) + y (017/1 +02772) - Vrj 3 , (1) 

where the Pauli matrices (rji) act in the sublattice 
(layer) space, and the Pauli matrices act in the val- 
ley space. Here, vp = y/Sta/2h is the Fermi velocity 
(vp ~ 10 8 cm/s), with t ~ 2.7 eV being the intralayer 
hopping matrix element. The k x (k z ) is the transver- 
sal (longitudinal) momentum calculated from the points 
K and K'. The potential difference between the layers 
opens up a gap 2\V\ in the bulk spectrum, while a spa- 
tial modulation, i.e. V — » V(x), breaks the translation 
invariance along the x-direction, thus only the total lon- 



E « t±. The effective Hamiltonian becomes 
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FIG. 2. The spectrum of the BLG structure for V /2 = 
50 meV and chirality = 0. The green area in the inset 
corresponds to the bulk spectrum. The mid-gap bound states 
for valleys K (full line) and K' (dashed line) have opposite ve- 
locities. The main figure shows the details of K-K' crossing 
region (shaded region in the inset). The curvature induced 
SOI shifts spin-up and spin-down levels in opposite directions 
by the SOI parameter f3. A magnetic field B assisted by inter- 
valley scattering A KK > results in the anti-crossing gap 2A g of 
two Kramers partners at k z = 0. If the chemical potential \i is 
tuned inside the gap [shaded area with ji « (Vb/2\/2)±/3], the 
system is equivalent to three Dirac cones (only one is shown in 
the main Figure), resulting in the helical mode regime. Here, 
P « 60 fieV (R = 5 nm), A KK , 30 /xeV, and A z = 30 /ieV, 
so the opened gap is 2A g « 30 /xeV~ 300 mK. 

The Hamiltonain Ho can be simplified for small volt- 
ages, \V\ <C by integrating out the A 2 and B\ degrees 
of freedom, which correspond to much higher energies 
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H = -V 7s ~ ( k l ~ k l) 7i - ^^k x k z T 372 , (2) 
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where the Pauli matrices 7^ act in the space of A\ and 
B 2 atoms. If the voltage changes sign at x = [for 
example, V(x) = —V(—x)\, this results in the closing and 
reopening of the gap. As a consequence, bound states, 
localized around x = 0, emerge within the bulk gap [15]. 
The eigenstates of Ho are characterized by k z and the 
valley degree of freedom r = zbl. For a step-like kink 
potential V(x) = (Vq/2) sgn(x) the energy spectrum is 
shown in the inset of Fig. 2. 

Now we include spin and aim at the realization of he- 
lical modes in BLG, which requires an analysis of the 
spin- full mid-gap states. At any fixed energy in the bulk 
gap, there are 2x4 states, where the factor 2 arises from 
spin-degeneracy. This means that the spectrum is topo- 
logically equivalent to four Dirac cones, each cone con- 
sisting of a pair of states with opposite momenta. On 
the other hand, helical modes are typical for systems 
with an odd number of Dirac cones. To effectively elimi- 
nate one Dirac cone at given chemical potential, the spin- 
degeneracy should be lifted by a magnetic field 5, giving 
rise to a new gap. Obviously, the opening of such a gap 
is possible only if there is level crossing in the system. 
The spectrum of the mid-gap states has support around 
K and K' . Therefore, if these points, projected onto the 
/c^-axis, are separated from each other, no crossing can 
occur. We thus see that the chiral angle is of a cru- 
cial importance for our purpose and the optimal choice 
is = (or very close to it). In this case, K 2 = YJ Z = 0, 
and the level crossing occurs in the center of the Bril- 
louin zone, at k z = 0, see inset of Fig. 2. We emphasize 
that in contrast to nanoribbons [16] the form of the edges 
of the BLG sheet does not matter provided the distance 
between edges and wire-axis is much larger than the lo- 
calization length £ of the bound state. 

Next, we allow also for spin-orbit interaction in our 
model. While the intrinsic SOI is known to be weak 
for graphene [11, 12], the strength of SOI in CNT is en- 
hanced by curvature [20, 23-25]. To take advantage of 
this enhancement, we consider a folded BLG which is 
analogous to a zigzag semi-CNT with = 0. All SOI 
terms that can be generated in second-order perturba- 
tion theory are listed in Table I of Refs. [20, 25]. From 
these terms only H so = f3r 3 s z is relevant for our problem; 
first, it is the largest term by magnitude, and second, it 
is the only term which acts directly in the A\-B 2 space. 
Here, Si is the Pauli matrix acting on the electron spin, 
and i = x,y,z. The value of the effective SOI strength /3 
depends on the curvature, defined by the radius R, and 
is given by f3 w 0.31 meV/R[nm] [20]. 

In the presence of SOI, the states can still be char- 
acterized by the momentum k Zl valley index r = zbl, 
and spin projection s = ±1 on the z-axis. The spec- 
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trum of Ho + H so can be obtained from the one of Hq 
by simply shifting E — > E — firs. This transformation 
goes through the calculation straightforwardly, and the 
spectrum in the presence of the SOI becomes 



E = /3ts± 



(hvpk z i 



± 



' (hvpk z 
4t± 



2V2. 



Vo_ 
V2' 



(3) 



The spin degeneracy is lifted by the SOI, giving a split- 
ting 2j3. As shown in Fig. 2, the level crossings occur 
between two Kramers partners at k z = 0: \K, t) crosses 
with I KM), and \K,l) crosses with \K',f). The KK'- 



crossing can occur provided \9\ < y 3(1 + y/2)t±Vo/^7rt. 
For the values from Fig. 2, we estimate this bound to be 
about 1°. As mentioned before, to open a gap at k z = 0, 
one needs first a magnetic field perpendicular to the SOI 
axis to mix the spin states, and second a K-K' scattering 
to mix the two valleys. Such valley scattering is described 
by the Hamiltonian H sc = A s kk ,ti + A^^/T^, where 

^kk' ) ^ s t ne scattering param- 
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eter for the bottom (top) layer of the BLG. The Zeeman 
Hamiltonian for a magnetic field B applied along the x'- 
direction is given by Hz = Azs xi with Az = g* /jLbB/2, 
where \±b the Bohr magneton. Here, g* is an effective g- 
factor due to the curvature of the fold and the localization 
of the bound state. Since s x * = s x cos ip+s y sin <p depends 
on x via the azimuthal angle (f(x) of the fold, we replace 
s x t by an average over the orbital part of the bound state 
wave function. This results in 2/tt < g* / g < 1, the 
precise value being dependent on the localization length, 
where g is the bare ^-factor of graphene. 

Using second order perturbation theory for j3 > A S KK , , 
Az, we find that the gap opened at k z = is given by 



*KK 



(4) 



see Fig. 2, which also contains numerical estimates for 
realistic parameters. We note that A^ is enhanced by 
electron-electron interactions [26], however, we neglect 
this supportive effect herein. If the chemical poten- 
tial is tuned inside the gap 2A^ [ji & (Vo/2\/2) ± /?], 
there are three right- and three left-propagating modes. 
Four states at finite momentum (two left-moving and 
two right-moving states) are only slightly affected by 
the magnetic field and thus can still be considered to 
carry opposite spins, meaning that the total spin trans- 
fer is close to zero and these modes are not contribut- 
ing to spin- filtering. In contrast to that, the two modes 
with k z « are helical modes: they have opposite ve- 
locities and opposite spins. Thus, similar to Rashba 
nanowires [18], the BLG quantum wire can be used as 
a spin filter device. 

Moreover, if the BLG is brought into proximity to an s- 
wave superconductor, the states with opposite momenta 
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FIG. 3. (a) The profile of the gate potential V{x) along the 
curved BLG. The reversed polarity at the two ends gives rise 
to mid-gap states, localized in x-direction. The density profile 
of three right-moving states, whose energy E is inside the 
gap A g , allows us to estimate the localization lengths £: (b) 
V 100 meV, E 29.13 meV, d = R (dashed white line 
in Fig. 2), f « 9 and 8 nm, and (c) V = 10 meV, E = 
3.49 meV, d = R, £ ~ 30 and 24 nm. We note that states with 
larger momenta have shorter localization lengths, (d) The 
localization length follows approximately £ = £ — £0 oc 1/y^Vo- 
The circles are extracted from our numerical calculations, £ = 
(x 2 ), for energies in the middle of the gap A g , equivalent to 
the white dashed line in Fig. 2. The lines are fits, £ oc 1/Vq P , 
for the states at large k z (p = 0.59), shown as dashed lines, 
and for the states near k z = (p = 0.52), shown as full lines. 



and spins get paired. Working in the linearized model of 
left-right movers [27], we obtain the effective Bogoliubov- 
de Gennes Hamiltonian for each of the three pairs, j = 
1, 2, 3, written in Nambu space, 



H 3 S = hvjkjXs + A s u; 2 X2, 



(5) 



where Vj is the velocity for the jth pair at the Fermi 
level and A s is the strength of the proximity-induced 
superconductivity, and the Pauli matrices \% (^i) ac ^ 
in the left-right mover (electron- hole) space. We note 
that we are in the regime corresponding to strong SOI 
where we keep only the slowest decaying contributions 
of the wave functions [27]. To determine the potential 
existence of MFs in the system, one can study the topo- 
logical class of Hi [28]. This Hamiltonian belongs to the 
topological class BDI. However, by analogy with multi- 
band nanowires [29] , additional scattering between states 
would bring the system into the D class. An alterna- 
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tive way of classification, which determines explicitly the 
number of MF bound states, is to study the null-space 
of the Wronskian associated with the Schrodinger equa- 
tion [30]. In our case, we find three MFs at each wire 
end in the topological phase defined by /S? g > + S/j 2 , 
where 5ji is the chemical potential counted from the mid- 
gap level A g . These MFs are generically hybridized into 
one MF and one non-zero energy fermion by perturba- 
tions such as electron-electron interactions and interband 
scattering. 

Numerical calculation. Above we have studied the 
system analytically, assuming a step-like potential. In 
this section we compare our results with the numeri- 
cal solution of the_ Schrodinger equation for the effec- 
tive Hamiltonian Ho + H so + H sc + Hz, with a more 
realistic (smooth) potential, V(x) = (Vb/2) tanh(x/d), 
where d is the distance between the gates. The spin- 
orbit interaction j3(x) is finite only within the curved 
region of the BLG sheet. Along the z-direction, the sys- 
tem is translationally invariant, so the envelope function 
is given by ty(x,z) = e lkzZ ipk z (x). The profile of ipk z (x) 
is presented in Fig. 3. The localization length follows a 
power law £ — £o oc 1/Vq 5 , with p « 1/2, and the shift 
£o < d is due to the finite distance between gates. In 
the limit d — » 0, where the analytical solution is ap- 
plicable, the localization length is essentially given by 
£ = 2 5 / 4 frup/V%£± [15] ? since corrections due to SOI 
are of negligible higher order in (3. 

Tunnel junction. The dependence of £ on the poten- 
tial Vo can be exploited to couple parallel wires. For 
instance, consider two similar quantum wires, running 
parallel to each other at a distance D. If £ <C D for each 
wire, then they are completely decoupled. However, low- 
ering the potential in both wires locally around a point 
zo on the z-axis, such that £o ~ D, we can enforce wave- 
function overlap, leading to a transverse tunnel junction 
between the two wires at Zq. In this way, an entire net- 
work of helical wires can be envisaged. We mention that 
such networks could provide a platform for implementing 
braiding schemes for MFs [31]. 

Conclusions. The confinement of states in BLG into 
an effective quantum wire is achieved by pairs of gates 
with opposite polarities, leading to eight propagating 
modes [15]. If the direction of the wire is chosen such 
that the chiral angle vanishes, both valleys K and K' are 
projected onto zero momentum k z . The SOI, substan- 
tially enhanced by curvature, defines a spin quantization 
axis and splits spin-up and spin-down states. A mag- 
netic field assisted by intervalley scattering opens up a 
gap at the center of the Brillouin zone. If the chemical 
potential is tuned inside the gap, three right- and three 
left-propagating modes emerge, so that the system pos- 
sesses helical modes, which are of potential use for spin- 
filtering. In the proximity to an s-wave superconductor, 
the BLG wire hosts Majorana fermions arising from the 
helical modes. By locally changing the confinement po- 



tential and thus the localization lengths, parallel wires 
can be tunnel coupled. This mechanism can be used to 
implement braiding of MFs in bilayer graphene. 

This work is supported by the Swiss NSF, NCCR 
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[1] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, 
M. I. Katsnelson, I. V. Grigorieva, S. V. Dubonos, and 
A. A. Firsov, Nature (London) 438, 197 (2005). 

[2] A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. 
Novoselov, and A. K. Geim, Rev. Mod. Phys. 81, 109 

(2009) . 

[3] R. Saito, G. Dresselhaus, and M. S. Dresselhaus, Phys- 
ical Properties of Carbon Nanotubes (Imperial College 
Press, 1998). 

[4] E. McCann, arXiv: 1205.4849. 

[5] A. K. Geim and K. S. Novoselov, Nature Materials 6, 183 
(2007). 

[6] R. T. Weitz, M. T. Allen, B. E. Feldman, J. Martin, and 
A. Yacoby, Science 330, 812 (2010). 

[7] J. R. Williams, T. Low, M. S. Lundstrom, and C. M. 
Marcus, Nature Nanotechnology 6, 222 (2011). 

[8] M. T. Allen, J. Martin, and A. Yacoby, Nature Commu- 
nications 3, 934 (2012). 

[9] A. S. M. Goossens, S. C. Driessen, T. A. Baart, 
K. Watanabe, T. Taniguchi, and L. M. Vandersypen, 
arXiv:1205.5825. 
[10] C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 226801 

(2005) . 

[11] H. Min, J. E. Hill, N. A. Sinitsyn, B. R. Sahu, L. Klein- 
man, and A. H. MacDonald, Phys. Rev. B 74, 165310 

(2006) . 

[12] M. Gmitra, S. Konschuh, C. Ertler, C. Ambrosch-Draxl, 

and J. Fabian, Phys. Rev. B 80, 235431 (2009). 
[13] S. Konschuh, M. Gmitra, D. Kochan, and J. Fabian, 

Phys. Rev. B 85, 115423 (2012). 
[14] F. Mireles and J. Schliemann, arXiv: 1203. 1094. 
[15] I. Martin, Y. M. Blanter, and A. F. Morpurgo, Phys. 

Rev. Lett. 100, 036804 (2008). 
[16] Z. Qiao, J. Jung, Q. Niu, and A. H. MacDonald, Nano 

Letters 11, 3453 (2011). 
[17] M. Zarenia, J. M. Pereira, G. A. Farias, and F. M. 

Peeters, Phys. Rev. B 84, 125451 (2011). 
[18] P. Stfeda and P. Seba, Phys. Rev. Lett. 90, 256601 

(2003). 

[19] M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 

(2010) . 

[20] J. Klinovaja, M. J. Schmidt, B. Braunecker, and D. Loss, 

Phys. Rev. Lett. 106, 156809 (2011). 
[21] J. Klinovaja, S. Gangadharaiah, and D. Loss, Phys. Rev. 

Lett. 108, 196804 (2012). 
[22] J. Alicea, Rep. Prog. Phys. 75, 076501 (2012). 
[23] F. Kuemmeth, S. Hani, D. C. Ralph, and P. L. McEuen, 

Nature 452, 448 (2008). 
[24] W. Izumida, K. Sato, and R. Saito, J. Phys. Soc. Jpn. 

78, 074707 (2009). 
[25] J. Klinovaja, M. J. Schmidt, B. Braunecker, and D. Loss, 

Phys. Rev. B 84, 085452 (2011). 
[26] B. Braunecker, G. I. Japaridze, J. Klinovaja, and 

D. Loss, Phys. Rev. B 82, 045127 (2010). 



5 



[27] J. Klinovaja and D. Loss, Phys. Rev. B 86, 085408 
(2012). 

[28] S. Ryu, A. P. Schnyder, A. Furusaki, and A. W. W. 
Ludwig, New Journal of Physics 12, 065010 (2010). 



[29] S. Tewari and J. D. Sau, arXiv:l 111. 6592. 
[30] J. Klinovaja, P. Stano, and D. Loss, arXiv:1207.7322. 
[31] J. Alicea, Y. Oreg, G. Refael, F. von Oppen, and M. P. A. 
Fisher, Nature Physics 7, 412 (2011). 



